Quasi-particle spin susceptibility in heavy-fermion superconductors : An NMR study 

compared with specific heat results 

Hideki ToiE 

Department of Quantum Matter, AdSM, Hiroshima University, Higashi-Hiroshima, 739-8530, Japan 

Kenji lshid£0 

Department of Physics, Graduate School of Science, Kyoto University, Kyoto, 606-8502, Japan 

Yoshio Kitaok^ 

Department of Materials Science and Technology, Graduate School of Engineering Science, 
Osaka University, Toyonaka, Osaka 560-8531, Japan 
(Dated: February 2, 2008) 

Quasi-particle spin susceptibility (x'^) for various heavy-fermion (HF) superconductors are dis- 
cussed on the basis of the experimental results of NMR Knight shift {K), NMR relaxation rate 
(l/Ti), and electronic specific heat (7e!) within the framework of the Fermi liquid model for a 
Kramers doublet crystal electric field (CEF) ground state. Xi is calculated from the enhanced 
Sommerfeld coefficient 'jei, and xti from the quasi-particle Korringa relation TiT{Kti)^ = const. 
via the relation of xti = {Na^J-b /Ahf)KTi where A^f is the hyperfine coupling constant, Na 
the Abogadoro's number and /is the Bohr magneton. For the even-parity (spin-singlet) su- 
perconductors CeCu2Si2, CeCoIns and UPd2Al3, the fractional decrease in the Knight shift, 
gj^obs ^ K°'"'{Tc) - K°'"'{T 0), below the superconducting transition temperature (T^) is due 
to the decrease of the spin susceptibility of heavy quasi-particle estimated consistently from x-y a-nd 
XTi . This result allows us to conclude that the heavy quasi-particles form the spin-singlet Cooper 
pairs in CeCu2Si2, CeCoIns and UPd2Al3. On the other hand, no reduction in the Knight shift is 
observed in UPts and UNi2Al3, nevertheless the estimated values of x-y and xti are large enough to 
be probed experimentally. The odd-parity superconductivity is therefore concluded in these com- 
pounds. The NMR Knight shift result provides a convincing way to classify the HF superconductors 
into either even- or odd-parity pairing. 

PACS numbers: PACS Number: 71. 28. -fa, 74.70.Tx, 75.40.Cx, 76.60.-k 



I. INTRODUCTION 

Extensive experimental and theoretical works for 
heavy-fermion (HF) systems have uncovered character- 
istic feature that a localized / state at high tempera- 
tures crosses over to a delocalized HF one at low temper- 
atures through the hybridization with conduction elec- 
trons. This occurs below a so-called coherent Kondo 
temperature Tk, being compatible with a renormalized 
heavy quasi-particle band width, leading to rich emergent 
phenomena such as either even- or odd-parity unconven- 
tional HF superconductivity (SC), anomalous magnetic 
or multiple ordering, etc. It has been confirmed that 
most physical quantities are described in terms of the 
Fermi-liquid theoryji e.g. revealing the enhanced Pauli 
magnetic susceptibility (x(T) ~ x(0)), the temperature 
(T)-linear coefficient in specific heat (C = JeiT) where 
"fei is the enhanced Sommerfeld coefficient, the T-square 
behavior in resistivity (p = po + ^T^), etc. 

As listed in Table I, the HF SC has been found in 
cerium (Ce) and uranium (U) based HF compounds, such 
as CeCuaSia,^ Celrlngfi, CeCoIn5^ UBcia,^''^ VPtsJ'^'- 
URuaSiajiS-iii^ UNiaAlg^iii^ UPdaAlaji^iiS-ii at ambi- 
ent pressure. In these heavy-fermion superconductors 



(HFS's), both large jei value and specific heat jump, 
AC/Cei{Tc) ~ 1, associated with the superconducting 
transition give unambiguous evidence that the heavy 
quasi-particle itself takes part in the formation of the 
Cooper pairs. 

In recent years, an intimate interplay between antifer- 
romagnetism (AFM) and SC has been the most inter- 
esting and outstanding issue in Ce-based HF systems. 
It is believed that the SC in CeCuaSia;^ Celrlns)^ and 
CeCoIngii emerges at the border of AFM even at am- 
bient pressure (P = 0). The finding of pressure- induced 
SC in CeCuzGesjia CePdsSiz,^^^ Celn3,2« and CeRhlngSi 
strongly suggest that AFM and SC are related to each 
other because P-induced SC occurs either when AFM 
vanishes or coexists with iti^ 

It is well known that Knight-shift measurements played 
vital role for establishing the Bardeen-Cooper-Schriefer 
(BCS) theory for spin-singlet s-wave superconductors pi 
and are the most powerful tool to identify which odd- or 
even-parity Cooper pairing state is realized in supercon- 
ductors. When the system undergoes a superconducting 
transition, the spin susceptibility Xs does decrease to zero 
at T=0 below Tc for even-parity (spin-singlet) supercon- 
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ductors as follows: 



= 2^ilNoY{T), 



where Y(T) is the Yosida function defined by^ 



Y{T) = - 



Jo 



(1) 



(2) 



and A'o is the density of states (DOS) at the Fermi level 
in the normal state, Nscsi^) is the DOS in the BCS 
superconducting state, and /(e) is the Fermi-Dirac func- 
tion. 

In transition metals and alloys, the electronic state is 
formed as follows. The crystal electric field (CEF) split- 
ting is much larger than the energy scale of the intra- 
atomic spin-orbit coupling (SOC), the orbital part of sus- 
ceptibility arises from the second order interband mix- 
ing effect with other bands, giving rise to T independent 
Van Vleck susceptibility. Accordingly, only its spin part 
yields the T dependence of susceptibility. Since both the 
Knight shift (K) and the susceptibility (x) depend on 
T, the contributions from the spin and the orbital parts 
are decomposed by taking the Clogston-Jaccarino plot of 
K(T) vs. x{T) with T as an implicit parameter:^ In the 
superconducting state, the spin part of the susceptibility 
decreases to zero in accordance with eqs. and 

On the other hand, the /-electron systems with the 
strong SOC are generally described by a total angular 
momentum, J ~ L + S . Furthermore, the CEF splits the 
lowest J manifold into several doublets and/or singlets. 
If the simple localized / electron picture can be applied 
for HF systems, the classical Van Vleck susceptibility be- 
tween the lowest lying CEF level forming quasi-particle 
bands and other CEF levels are given by 



Xvv = 2Na 
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where Na the Abogadoro's number, gj the Lande fac- 
tor, fiB the Bohr magneton, Z the distribution function. 
Ea and E/^ are ground (a) and excited (/3) state ener- 
gies. Note that this interband contribution reveals a T 
dependence when the CEF splitting is compatible with 
T. That is, the T dependence of both K and x are not 
due to real spin, but to fictitious spin including the inter- 
band contributions. Therefore, it is difficult to conclude 
the odd-parity superconducting state in case of the Van 
Vleck contribution being dominant )2242iii even though 
K unchanges across Tc- Anyway, in the HF systems, 
it is not so clear to extract the quasi-particle spin sus- 
ceptibility via the Clogston-Jaccarino plot, and it is 
crucially important to reconfirm the x"^^ due to heavy 
quasi-particles. 

In this paper, we show that the quasi-particle spin sus- 
ceptibility x"^^ in HFS's is reasonably estimated from the 
NMR and specific heat results on the basis of the Fermi- 
liquid theory with a Kramers doublet CEF ground statei 
In § 2, we remark that x''^ is independently estimated 



from the Sommerfeld coefficient and the NMR relax- 
ation rate 1 /Ti , which are consistent with estimated 
from the Clogston-Jaccarino plot. In § 3, we describe 
the basic assumptions for the present analysis. Section 
4 is devoted to the analysis of the experimental data for 
various HFS's in terms of the Fermi-liquid model. We 
also discuss the validity of the present model. After the 
brief comments on the quasi-particle susceptibility in the 
superconducting state in § 5, we discuss the relation be- 
tween the Knight shift and the parity of the order pa- 
rameter (OP) in the superconducting state in § 6. From 
these analyses, CeCu2Si2, UPd2Al3, and CeCoIns are re- 
inforced to be an even-parity (spin-singlet) superconduc- 
tor, whereas UPta and UNi2Al3 are reconfirmed to be an 
odd-parity (spin-triplet) superconductor. Finally, § 7 is 
devoted to several concluding remarks. 



II. NMR PARAMETERS AND THEIR 
RELATION TO FERMI-LIQUID PARAMETERS 

A. Quasi-particle susceptibility 

In the Fermi-liquid theory for HF system, the elec- 
tronic specific-heat coefficient jei (in units per mole) and 
the quasi-particle spin susceptibility x"^^ (in units per 
mole) are enhanced by the mass enhancement factor 7, 
and the magnetic enhancement factor xq associated with 
the contribution of /-electrons to the spin susceptibility, 
respectively. These are approximately written asm 



X 
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NA9jJeffl^BPH£F)X0, 



(4) 
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where fc^ is the Boltzmann factor, Je// the effective spin, 
{ef) the one-spin DOS at the Fermi level for bare /- 
electrons. Here, the DOS per spin of quasi-particles is 
approximately written by 



From eqs. Q) — ||HJ), x'^^ is related to jei as 



(6) 



Xo 



NA9jJ^ff^i'BP*(.SF)^ 



leigPtffPB 



r2 1,2 



(7) 



where R is the so-called Wilson Ratio, defined by i? = 

When an external magnetic field is applied, x'^^ Pro- 
duces a small additional magnetic field at the nucleus. 
This gives rise to a Knight shift {K'^p) due to the quasi- 
particles that is scaled to through the hyperfine- 
coupling constant A^f (in units per /zb) at q = as 
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TABLE I: Physical quantities and NMR parameters for various HFS's. Tk, "fei, Tc, A^^f, and (l/Tir)Tc,i are Coherent Kondo 
temperature, Sommerfeld coefficient (at Tc), superconducting transition temperature, hyperfine couphng constant, NMR spin- 
lattice relaxation rate divided by temperature at Tc, respectively. Here i =|1,-L (see text). 
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From eqs. Q and ijHJl, we obtain the quasi-particle spin 
contribution to the Knight shift, K'^p, as 



KIP = 



R. 



(9) 



Here it is noteworthy that no interband contributions are 
included in this formula. 



B. Quasi-particle Korringa relation for isotropic 
system 

The nuclear spin-lattice relaxation process basically 
occurs through the spin-flip process between the state 
I — 1/2} and | -j- 1/2) of conduction electrons. Provided 
that HF state is realized at low temperatures below Tk, 
K''P is also extracted using the NMR relaxation rate, 
l/Ti. 

In general, l/Ti is given bji^ 

where 77V is the nuclear gyromagnetic ratio, A^y(q) 
the q-dependent hyperfine-coupling constant and 
Imx± (q, the imaginary part of the transverse 
component of dynamical susceptibility x(q, w) of 
quasi-particles. If the HF state is formed above Tc, 
the imaginary part of x(Qi is expressed within the 
random phase approximation (RPA) as^*^ 



(11) 



where XQ(q, w) is the dynamical susceptibility of /- 
electrons and x{ci, w) an enhancement factor due to the 
electron-electron interaction as 



x(q,^) 



(12) 



Since the imaginary part of Xo (q^ near at the Fermi 
level ef is expressed asm^^ 



(13) 

then we have 



q \ q / 

(14) 

where Xq = x(qj'^ ~^ 0) at the NMR frequency and the 
bracket represents averaging over the Fermi surface. If we 
are neglect the q dependence of the hyperfine coupling, 
then Ah / (q) is replaced by Ah / and the NMR relaxation 
rate 1/Ti is expressed as 



1- = ^'L^^^nAhs)\gjJcfffpf{eFfkBTl^xVj , 

(15) 

From eqs. (jSJ, JHl and (|15|l . we have 
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where IC{a) = \ J2qXq) /^o is the enhancement fac- 
tor and a is the parameter associated with the Storner 
factor If the enhancement factor (j2q Xq^ is the same 

as that at q = 0, i.e, JC{a) = 1, this relation is similar to 
the so-called Korringa law. Then, when ferromagnetic 
and antiferromagnetic spin fluctuations become domi- 
nant due to electron-electron correlation effect, K.{a) < 1 
and IC{a) > 1 are expected, respectively. Here we should 
notice that 1 /Ti does not include any interband contribu- 
tions associated with the excitations across the CEF split- 
ting. This is because 1 /Ti involves the imaginary part of 
the dynamical susceptibility in the low energy limit as 
w — > 0, as seen from eqs. (|10ll and H14|) . Anyway, K'^p 
derived from the quasi-particle Korringa relation does in- 
deed give the quasi-particle susceptibility. 



III. BASIC ASSUMPTIONS FOR ANALYSIS 

A. CEF ground state 

In Table I, we notice that the NMR l/Ti for HFS's 
is enhanced by more than ten times in comparison with 
that for non-HF materials,^-'^i^5i^7i^8,29,33,34,35,56,57,58 

instance 1/TiT = 33 (1/sK) for ^^^Pt of pure Pt metalS, 
and 1/TiT = 0.033 (1/sK) for ^"^Al for LaPd2Al3,^ etc. 
This means that the heavy quasi-particles are involved 
in the NMR relaxation process. Since the non-Kramers 
ground state is magnetically intact, it is difficult to ex- 
plain the enhanced relaxation rates observed in HFS's by 
the present Fermi-liquid model. As mentioned in the pre- 
vious section, the spin-flip process of conduction electrons 
would not occur through the Van Vleck type inter-band 
excitations in the low-energy limit as w — > 0. Thus, we as- 
sume that the enhanced relaxation rate is ascribed to the 
magnetically degenerated HF ground state. Namely, we 
assume that several doublets and singlets lying in lower 
energy region than I^bTk are renormalized into a ground 
state having Kramers degeneracy via the c — f mixing 
effect. 



B. Effective moment below Tk 

In the spherical system without any CEF, the effective 
spin is defined as J^^-j = J( J + 1) and the effective mo- 
ment is ^eff = gj\/ J{J + ^}fJ'B- For the free electrons, 
J^ff = y/3j2 and fieff = 1.73/iB. In HFS's, the CEF 
level scheme, that is determined experimentally, should 
be taken into account. Furthermore, several doublets and 
singlets are expected to be renormalized by the c— / mix- 
ing effect in the heavy-Fermi-liquid state. It is, however, 
underlying issue to estimate to what extent the CEF level 
splitting affects x"*^ in HFS's. 

For Ce-based HFS's having the clear CEF level split- 
ting, the basis function of a CEF ground state is not 



spherical but anisotropic. In order to introduce the 
anisotropy, we assume that the effective moment at low 
T's is attributed to only the CEF ground state level and 
we do not take the average of each direction. Then, 
the effective moment is tentatively replaced by fJ-l,ff = 
gjJlffHB = gj{0\Ji\0)^iB where |0) is the basis func- 
tion of a CEF ground state. For CeCu2Si2, using the 
CEF scheme, the respective effective moments at low 
T for _ff ± c and H \\ c are estimated as = 

9j1^b{0\J±\0) ~ 0.90Ais and = gjtiB{0\J\\\0) ~ 

1.07/iB. For CeCoIn5j^^ the effective moment for ± c 
is also calculated as fi-^j^ ~ 1.28^b. 

By contrast, for U-based systems, a clear CEF level 
splitting has not been reported yet. This has been gen- 
erally accepted as the itinerant nature of U 5f-electrons. 
So that, the basis function would be rather spherical. In 
our analysis, the effective moment for U-based systems 
is tentatively taken as /ie// — 1.73/is, i.e., free electron 
value. The magnetic anisotropy is, however, expected 
from the various magnetic measurements. Therefore the 
effective moment for the magnetically easy axis is taken 
to be Me//' ~ 1-73/iB, while that for the magnetically 

hard axis is to be n'^f}'^ = 1.73 ^^B x ^ xtarJ xfasy where 

the ratio, Xhard/Xeasy^ is estimated from the anisotropy 
in magnetic susceptibility at listed in Table I. 



C. Hyperfine coupling constant 

The q-dependent hyperfine couphng constant Ahf{q) 
in eq. H10|l is the spatial Fourier transform of the hy- 
perfine coupling Ah fir) between nucleus and Ce/U ions 
separated a distance r. If the nucleus is coupled to 
only one Ce/U ion, Ahf{q) is independent of q. In 
the HFS's, the calculation of A/j/(q) would require one 
to sum over contributions from the next nearest Ce/U 
neighbors. However, Ahfir) is not so simple. Therefore, 
under the most naive assumption, we neglect the wave 
number dependence of the hyperfine coupling constant, 
i.e., Ahfiq) = Ahj. 



D. Anisotropic quasi-particle Korringa relation 

In HF systems, the magnetic anisotropy is inevitable 
even in the HF state because the HF state at low T is 
formed through the c — f hybridization. It should be 
noted that a large magnetic anisotropy in xl^ii — c- 
axes) arises from the anisotropy in the hyperfine-coupling 
constant the effective moment (iJ'lff) and the en- 

hancement factor (xo)- Unfortunately, no anisotropic 
Fermi liquid theory have been provided yet. Further- 
more, the CEF splitting and the c — f mixing effect 
are not so clear in quantitative level. Therefore, from 
a phenomenological point of view, a uniaxial magnetic 
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anisotropy is taken into consideration as 



(17) 
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NaIJ.1 
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where i 



and _L which mean the respective parallel 



and perpendicular components to the field direction, e.g., 
X|| = Xc and x± = Xa = Xb- Then the Knight shift in 
eq. (jSJl or (O is replaced by anisotropic Knight shift as 



k: 



qp _ 
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hf qp 



The quasi-particle Korringa relation, eq. (fT?))l . is modi- 
fied using anisotropic parameters, ^Jjy, Me//i and ICi{a), 
as 



1 \ Snks 1 1N^\ fr^qp-.2r' i \ 
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IV. DISCUSSION : THE NORMAL STATE 
A. Quasi-particle Korringa relation 

With above assumptions, we calculated the effective 
values of the quasi-particle Knight shifts at Tc from 
eqs. (fTH|) — by using parameters 7e;, ^hf^ ^^'^ 

{l/TiT)T^,i listed in Table I. Here we defined the effective 
quasi-particle Knight shift associated with ^ei including 
the Wilson ratio R as K^^i = K^^ /Ri while that associ- 
ated with Ti including the enhancement factor K,i{a) as 
Kr^^i = Kl^^JlCi{a). Thus calculated values K^^i and 
Kxi^i for various HFS's are listed in Table II. 

As clearly seen in the eqs. (|18ll = (|20|l . both \K^^i\ and 
T^Ti.il depend on the size of /ie//- In order to check the 
vahdity of the present analysis, |-ft'Ti,i| (= \Kf\\/K.i{a)) 
is plotted against \K^^i\ (= \Kl^\/ Rt^) for various HFS's 
in Fig. 1. Most of valu es are consistent with a relation 
of \KTi^i/ K^^i\ = Rt^ y/ ICi {a) = 1 except CeCu2Si2 and 
CeCoIns. 

According to the Fermi-liquid theoryji^ the Wilson 
ratio varies in the range of li?2. Simply, i? = 1 is the 
limit of electron-electron correlation parameter J7 — s- 0, 
while i? = 2 is the limit of [/ oo. The recent the- 
oretical studies have predicted that the mass enhance- 
ment factor is comparable to the magnetic one, i.e., the 
Wilson ratio is i? « 1.2i2ii^ Actually in CeCuaSia, un- 
der the present assumption of fj,eff{~ ^^J■B), we obtain 
Rt^ = 1.24 from eq. {Tj) with the experimental values, 
-fei" ^ 1000 mJ/(moleK2),2 and x°'"iTc) w 6 x 10"^ 
(emu/mole) A good agreement between the theoret- 



ical and the experimental values strongly suggests the 
Wilson ratio in the HF liquid state is on the order of 
unity. 

For U-based HFS's, on the other hand, the absence of 
the clear CEF splitting prevents quantitative estimates 
for /ie//) consequently for Rt^. However Rt^ would be 
more close to unity than that in Ce-based systems be- 
cause U-based HFS's share more itinerant nature than 
Ce-based ones do. Thus the relation of \KTi.i/Kj^i\ — 
Rt^ \/ ICi (a) = 1 suggests that /^e// ~ 1.73/iB is proba- 
bly a good approximation for the effective moment in the 
U-based systems. Instead, if we assume /^e// = l/^s and 
Rt^ = 1 in U-based systems, we obtain /Ci(a) « 4 that is 
comparable to that in Ce-based HFS's in which antifer- 
romagnetic spin fluctuations are dominant and systems 
are close to a magnetic quantum critical point. ^^■^■^ This 
is, however, unlikely to occur in U-based HFS's becasue 
l/TiT = const, behavior strongly suggests that antifer- 
romagnetic spin fluctuations are not dominant at low T's. 

To be emphasized here is that the quasi-particle Kor- 
ringa relation is actually hold in the HFS's, though it is 
still unclear to evaluate R^^ and /Ci(a) experimentally. 
We also calculated the susceptibilities X7,i(= xT l^i) 
and XTt,i{= xT\/^i{<^)) at Tc by using eq. CHI) (see 
Table III). 



B. Van Vleck susceptibility 

Here we deal with the Van Vleck susceptibility Xvv,i- 
According to the theoretical argumentsr'^Sii^ii^ Xvv.i is 
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TABLE II: Observed reduction of the Knight shift 5K°^' 
K^p/Rt^ and Kti 



K°'"'{Tc) ~ K°'"'{T ^ 0) at T -> K for various heavy fermion 



superconductors. Kj 



K'lP^K.ia) are the calculated effective Knight shifts in the normal state from 
the electronic specific heat (7ei) and the NMR relaxation rate (1/Ti) assuming Ahf{q) = Ahf, respectively. Here we assume 
fj,eff ~ 1.73//S for U-based HFS's (see § 3). For UPts, the respective 5K°^'"s for C and B-phases are listed. 
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and XTi,i provide the quasi-particle susceptibility jR 
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the quasi-particle Korringa relation {TiT{K'''')^ JC{a) — 



cor.stant) and the electronic specific heat (if''^ oc 7e;ii') just 

for 
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ab ove Tc, respectively. — Here we assume Ahf{q) — Ah f 
calculations of Kt^, and /Xe// ~ 1.73/xs for U-based HFS's 
(for details, see § 3). The solid line shows R^JlCia) = 1. 



and xT \/ ■ If assume Rt^ « 1, x' 
That is, the Van Vleck susceptibility is approximately 
written as Xvv,i — x'i"^{Tc) — X-y,!^ The calculated val- 
ues of xvv,i are listed in Table III. Note that Xi(T'c)'s 
are larger than x-y,i and XTi,i except for UBeiajS^ and 
Xvv,iS are comparable with or larger than the quasi- 
particle susceptibility. This seems to be consistent with 
the theoretical prediction that xvv can become compa- 
rable with or larger than the quasi-particle susceptibility 
(X«P) via c - / niixing,2MMi42ii4^ This in turn indi- 
cates the observed Knight shift K°'^^ involves the Van 
Vleck contribution. 



QUASI-PARTICLE SUSCEPTIBILITY IN 
THE SUPERCONDUCTING STATE 



Here, we discuss about a relation between the Knight 
shift and the parity of the superconducting order param- 
eter (OP). In the superconducting state, the DOS for 
quasi-particle excitations is written asS 



P*{eF)\e\ 



generally expressed as, 



XVV.i = XT,i 



qp 

At 1 



(21) 



jO*(£f ) |£| 
47r 



27r /•TT 



sm 



(22) 



where xta is the total susceptibiUty and xta = Xi (Tc) 
at Tc- As discussed in the previous section, both X'y,i 



where the integral is over the solid angle of 17k and p*{£f) 
is the quasi-particle DOS at the Fermi surface in the nor- 
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mal state defined by eq. 
tion and is given by 



©. Here A(k) is the gap func- 



A(k) oc 



l*(k)| 



(/ = even), 
d(k)|^ ± |d(k) X d*(k)| (/ = odd), 



(23) 

where ^'(k) is a single even function for a spin-singlet 



pairing state (5 = 0,/ =even), whereas d 
a d- vector for a spin-triplet pairing state {S — 1, 1 —odi 
Due to the k dependence of the OP, the superconducti 
energy gap vanishes at points and/or along lines on t 
Fermi surface. 

The quasi-particle spin susceptibility, x'^^i below Tr 1 
even-parity superconductors is expressed as 



and 



Y{T) = - 



(2 



(2 



Thus decreases down to zero as T ^ for the eve 
parity (spin-singlet) superconductors, when the surfa( 
and/or impurity effect are absent. Notice here th 
the anisotropy of the effective moment A'e// causes 
anisotropic decrease of x'^^ below T^,, depending on a fie 
direction. 

Differently from an even-parity pairing state, 
an odd-parity pairing state, the Zeeman interactic 
(l/2)x''^|H - dp, due to the spin of the Cooper pa 
(d-vector) should be taken into consideration. Unle 
d-vector is locked to the lattice, the application of fie 
makes the d-vector rotate as d_LH so as to minimize t 
Zeeman energy. In this case the quantization axis for t 
spin of the Cooper pairs is always parallel to the ma 
netic field. Namely, the quasi-particle susceptibility \ 
is essentially the same as that in the normal state; x 
does not change below Tc. If the spin of the Coop 
pairs is "frozen" with the lattice^SiS it is expected th 
the quasi-particle susceptibility, x'''(7^)j decreases wh 
H II d even for an odd-parity pairing state. It has be 
argued that the strong intra-atomic SOC may make 
vector freeze. But, unfortunately, the inter-atomic SC 
for the pairs remains quite obscurei^ 



VI. DISCUSSION: THE SUPERCONDUCTING 
STATE 

A. Knight shift and parity of the Cooper pairs 

Figure 2 shows the T dependence of the Knight shift 
for CeCuzSia^ CeCoIna^ and UPdzAlaiSLSa Figure 3 
shows the T dependence of the Knight shift for UPta 
(for the C-phase in multiple phases)^ii2S, UNi2Al3'^^, and 
Sr2Ru04.''^ As clearly seen in Fig. 2, the reduction of the 
Knight shift below Tc defined by 



is independent of the crystal direction for CeCu2Si2, 
CeCoIns, and UPd2Al3, ^3,25,27,28 ^he fractional decrease 
in K°^^ below Tc is ascribed to the reduction of due 
to the formation of the spin-singlet Cooper pairs. In Fig. 
3, on the other hand, the 5K°^'^^s unchange across Tc for 
UPta and UNi2Al3, and also for Sr2Ru04> 



31.32.33 




T/T, 




T/T„ 




T/T„ 



Fig.2 

FIG. 2: (a) T dependence of the "^CM-Rwi^itlsMfesYfdPtSf^ia 
ented powder sample measured by Ueda et a?.(R G f K^i ). (b) T — 



5K' 



obs 



K°'"{Tc) - K°'"'{T 0), 



(26) 



dependence of the ^'^AI-Knight shifts measured by Kyogaku 
et al. (Refsi2£i2^). Closed triangle shows the Knight shift 
for _ff _L c in oriented powder sample. Closed and open circles 
show the respective Knight shifts for _ff _L c and || c in sin- 
gle crystalline sample, (c) T dependence of the ^^^In-Knight 
shifts for oriented powder sample measured by Kohori et al. 
(Reff— ). The solid lines in figures are calculations using the 
Sommerfeld coefficient y^i for a d-wave singlet Eig model (see 
§ 6.1). Here we assume /Xe// ~ 1.73/iB for U-based HFS's (for 
details, see § 3). 



The solid lines in the figures are calculated from 
eqs. ©, and ((23), with the anisotropic energy 

gap A(T, k). In our calculations, we tentatively choose 
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the gap function as d-wave Eig 



52.67 



A(T,k) = A{T)ik,±iky)k, 
= A(r)sin6'cos6'e^*'^ 



(27 



with the T dependence of the gap which is the same a 
the BCS theory as 



A(r) = A(0)tanh 



vrfesT, AC /T, ^ 



, (28 



where AC /C is the specific heat jump, and a — 
and A(0) is the magnitude of the superconducting en 
ergy gap. Since the energy gap for the Eig symme 
try vanishes along the hne on the Fermi surface, i.e 
pl{e) oc e {e ^ A(0)), it can explain the behavio 
of NMR relaxation rate observed in the HFS's, wher 
1/Ti oc /o°°p*(£)V(£)(l - f{e))de oc for line-nod 
gapped superconductors. Furthermore, we use 
K'^P assuming Rt^ — 1. This is because Kt^ is some 
times enhanced by IC{a) for Rt^ — 1, so that the spr 
part may be overestimated. 

The calculations reproduce the observed fractional de 



obs 



0.12 % and SKf'' < -0 



crease in K , i.i 

% for CeCuaSis (^^Cu NMR)^^; SK"'" « 0.08 - 0.1 
% for UPdaAla (^^Al NMR)2LS; SK"'"' « 0.5 - 0.6 9 
for CeCoIng ("^in NMR)^^ Here, the T dependence c 
calculated Knight shift was obtained by using the ap 
propriate energy gaps of 2A{0)/kBTc = 5, 4, and 7 fo 
CeCu2Si2, UPd2Al3, and CeCoIng, respectively. Any 
way, it is demonstrated that x'^^'s decrease to zero a 
T — > due to the formation of spin-singlet Cooper pair 
in CeCu2Si2, CeCoIng, and UPd2Al3. As a matter c 
fact, x-y (also xti) provides a reliable way to deduce x' 
in the HF state. 

Likewise, in Fig. 3, the and Kt^ in UPta arl 
estimated to be very large, which are compatible with 
the previous estimates both by the Clogston-Jaccarino 
plot and the high-T Curie- Weiss fitting to the Knight 
shift.— Also those in UNi2Al3 are comparable to the val- 
ues in UPd2Al3. Similar result is obtained for Sr2Ru04 in 
which a calculated value of x-y (^7) is in good agreement 
with the spin part of the susceptibility estimated from 
the Clogston-Jaccarino plot reported previously (See Ta- 
bles II and III)i2i No significant reduction in K°'"^ is 
observed in these compounds, ruling out a possibility 
for even-parity (spin-singlet) state. The solid lines in 
the figures are expected if a spin-singlet pairing states 
were realized with 2A{0) /ksTc =10, 5, and 4 for UPts,^'' 
UNi2Al3)2^ and Sr2Ru04r^ respectively. These experi- 
mental results reinforce that UPta, UNi2Al3, Sr2Ru04 
are the odd-parity superconductors, where the d-vector 
rotates freely as d _L H in the C phase in UPta;^ while 
it does as d II c-axis (the odd-parity state with paral- 
lel spin pairing lying in the ab plane) in UNi2Al3^ and 
Sr2Ru04.^1 

To make discussion more clear, we define the normal- 



-15 
-10 



-5 - 



1 . . . , 1 1 . 
(a) UPt, 


. 1 1 1 1 . 
H//a (•), 




H//b(A) _ 




/ H//C 


OCDCDO CDD CrOOC 








1 1 1 1 1 1 1 1 


, 1 . . . , 



0.5 



T/T^ 




T/T, 



T/T. 



FIG. 3: (a) T dependence of the H^feiCiKgto^tftsffcKiYSSka 
ious field direction (Refi— ). (b) T dependence of the ^^Al- 
Knight shift for H J. e-axis (Ref ,^ ). (c) T dependence of the — 



^''Ru- Knight shift for H _L c-axis (Refp^). The solid lines in 
both figures are calculations using the Sommerfeld coefficient 
7e! for a d-wave singlet Eig model (see § 6.1). Here we assume 
fieff ~ 1.73/iB for U-based HFS's (for details, see § 3). 



ized residual susceptibility at T ^ as 



S 



SK' 



obs 



(29) 



where we assume 5K°^^ — {A^f /Naiib)5x°'^^ ■ In Fig- 
4, the residual susceptibility ratio 5{%)= {x-j{Tc) — 
Sx°'"^)/x-f{Tc) is plotted against Xj- Note that S is al- 
most zero (or below zero) for CeCu2Si2, CeCoIng, and 
UPd2Al3, ensuring that the x"^^ in these compounds de- 
creases down to nearly zero at T — > regardless of the 
crystal direction. Together with the fact that 1/Ti oc 
below Tc, which excludes the s-wave BCS state, a d- 
wave type anisotropic superconducting state is realized 
in CeCu2Si2, CeCoIng, and UPd2Al3. 

By contrast, S « 100 (%) for UPt3 and UNi2Al3 gives 
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FIG 

as S 
is re.ated 



150r 



St.RuO, UNi,Al,(±) UP1,(1) UBe,-(av) 

100 A ® A — m ^ 9 ^ 

; URij,Si,(±l URu,St( II } UP1,( II ) 



50- 



_CcCu,Si,(||) 
^ al 0.5T, 



UP(iAl,(l)n QCeCu,Si,(Sil) 
OceCaIn,(±) 



10' 10' ^ 10° 10' 

X , (X 10' emu/mol) 



4: Normalized residual susceptibility at T 
= {X-,{T.) - Sx^'l/XiiT.) plotted against x-,{Tc 

to the quasi-particle pj^^ieptibility as x''^(Tc 
. For UPts, data points HrgFeupKlclBtlflibed'. fcha6kie val- 
ues ifi B-phase. 



d|elined 
which 



evidence for an odd-parity (spin-triplet) superconducting 
state. It was reported that K"'"' unchanges across in 
UBeia^ and URu2Si2 as welli^ Although an odd-parity 
pairing state is most likely for URu2Si2 and UBeia from 
the present analyses, further precise measurements for 
high-quality single crystalline samples are desirable. Also 
the present analysis is applied to Sr2Ru04, reinforcing 
that it is the spin-triplet superconductor 



B. Comments on quasi-particle susceptibility in 

UPt3 

Here we discuss the quasi-particle susceptibility in 
UPta. In the previous paper we reported that the 
value of Sxb'"' ~ ^Xc'"' « 1 x 10~^ emu/mole is two orders 
of magnitude smaller than the value of measured suscep- 
tibility which is comparable with the present estimates, 
Xj and xti- Fig. 5 shows the T dependence of the ^^^Pt 
Knight shift {^^^K) in UPts at various magnetic fieldsiS 
As seen in the figure, when the field strength is decreased, 
the Kf" for the c-axis (Fig. 5(c)) and K°''' for the b- 
axis (Fig. 5(b)) decrease across Tc aX H < b kOe and 
H < 2.3 kOe, respectively. By contrast, Kf"" (Fig. 5(a)) 
unchanges across Tc down to the lowest field oi H = 1.764 
kOe. Remarkably, ^^^K depends on the strength of mag- 
netic field and its crystal orientation. These results can- 
not be explained by the even-parity (spin-singlet) pairing 
model at all, but is consistent with the odd-parity (spin- 
triplet) oneM 

In order to discuss the quasi-particle susceptibility in 
UPta within the present model, we focus again on the 
NMR relaxation behavior. It should be noted that 1/Ti 
is enhanced in the HF state via the transverse component 
of dynamical magnetic susceptibility which is also respon- 
sible for the mass enhancement ^2*^ In most uranium 
compounds except UPd2Al3,^^ it is widely believed that 
a 5/^ configuration may lead to a non-Kramers ground 



state. For the HF state involving only a singlet CEF 
ground state, however, l/Ti may not be enhanced, as 
mentioned in § 3-1. The fact that l/Ti is enhanced in 
HFS's suggests that the Kramers doublet or magnetic 
CEF ground state is responsible for the formation of HF 
state, associated with more itinerant nature of uranium 
5/ electrons than of cerium 4/ electrons. In this context, 
the reason why the reduction of 6x°''^ in the SC state 
is so small may be ascribed to the incomplete lock of d- 
vector to the crystal, that is, the spin degree of freedom 
of the Cooper pairs is not always completely locked to 
the crystal axis. This suggests that the Knight shifts for 
Kj^''^ and K°''^ decrease markedly at further lower fields. 
The ^^^Knight shift measurements at fields lower than 
1.7 kOe are needed to check thisi^ 



(a) UPtj #3s H//a 



(b) UPtj #3s H//b 



(c) UPtj #4 H//c 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 
T (K) T (K) T (K) 

Fig.5, 

H. Tou, K. Ishida, Y. Kitaoka 

FIG. 5: Temperature dependence of the ^^^Knight shift in 
UPts at various magnetic fields for {a,)H \\ a, {h)H \\ b, {c)H || 
c (RefiS). Arrows show Td and Tc2- 



Meanwhile, recent theoretical studies proposed that 
the origin of the tiny reduction of 6x'^^ is possible to 
be explained by assuming that the U(5/^)-derived non- 
Kramers singlet CEF ground state is hybridized with 
conduction electrons>£^i£^ It was argued that the x'^ for 
a singlet CEF ground state is not enhanced, whereas 
the DOS of quasi-particles is enhanced via the on-site 
correlation effect giving rise to the heavy-Fermi liquid 
behavior i24 It was predicted qualitatively that the NMR 
1/Ti is enhanced by an order of 1/T|-, while x'^^ is not 
enhanced by correlation effect. Unfortunately no direct 
evidence for the CEF splitting in uranium compounds 
prevents the complete understanding of the HF state for 
the 5/^ configuration. Further quantitative discussions 
in theory in addition to experimental efforts are needed 
to test whether or not this theoretical model is applica- 
ble to UPta as a realistic model. Anyway, the physical 
properties of the heavy-Fermion having a singlet CEF 
ground state should be addressed by both experimental 
and theoretical works in future. 
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TABLE III: Reduction of the susceptibility ^x"*" = {Na^b /Ahf)5K°^' at T ^ OK for various heavy fermion superconductors. 
The calculated effective susceptibility, x-i — /Rt^ and Xti = x''^ \/^{'^)^ by using the relation of X7,Ti ~ {NaI-^b /Ahf)Kj^Ti 
assuming Ahf{q) = Ahf- Here we assume /tie// ~ 1.73/is for U-based HFS's (see § 3). For UPts, the respective 5x°'""s for C 
and B-phases are listed. The Van Vleck susceptibility is calculated from xvv = X°'"'{Tc) — x-t- Xc's are the spin suceptibilities 
for pure metals. 





nucleus orientation 


Sxf' 




XTi.i 


Xvv,i 


Xc 


Parity 








(10"2emu/mol) (10" 


■^emu/mol) (10" 


"^emu/mol) (10" 


"^emu/mol) (10" 


"^emu/mol) 




CeCu2Si2 


^^Cu 


II 


>0.26 


0.52 


1.88 


1.43 


0.001 


even 




s^Cu 


± 


0.44 


0.36 


1.35 


1.21 








2«Si 


II 


— 


0.52 


1.6 


1.43 


— 






2^Si 


± 


0.37 


0.36 


0.36 


1.21 






UBei3 




Average 





1.50 


1.89 


0.0 


0.002 


odd(?) 


UPts 


195 


II 


0/0.0126 


0.31 


0.23 


0.14 


0.009 


odd 




196pt 


± 


0/0.0099 


0.57 


0.67 


0.28 






URu2Si2 


29Si 


II 





0.089 


0.114 


0.371 




odd(?) 




2^Si 


± 





0.029 


0.025 


0.121 






UPd2Al3 


27 Al 


II 


0.13 


0.11 


0.095 


0.15 


0.002 


even 




27 Al 


± 


0.204 


0.205 


0.177 


0.65 






UNi2Al3 


2^A1 


II 




0.07 


0.17 


0.118 


0.002 


odd 




27 Al 


± 





0.17 


0.17 


0.34 






CeCoIns 


"5in 


II 












even 




"^In 




0.32 


0.27 


0.94 


0.43 






Sr2Ru04 









0.0535 


0.068 


0.0415 




odd 



VII. CONCLUDING REMARKS 

Wc have presented how the spin susceptibility x''' 
of quasi-particles in the HFS's is reasonably estimated 
within the framework of the Fermi-liquid theory, focus- 
ing on the NMR and spccific-hcat results. The present 
analysis for the Kramers doublet CEF ground state has 
provided two important results on the HF superconduc- 
tivity: 

1. In the normal state, using the value of TiT= const, 
and the T-linear coefficient 7e; in specific heat, the 
respective values, Kt^ and K^, are related to the 
quasi-particle Knight shift K'^'p . The quasi-particle 
Korringa relation was found to be valid far below 
Tk in the HFS's. For U-based HF superconductors, 
it was expected that any significant wave-number 
dependence would be absent in low-lying excita- 
tions. This means that the on-sight correlation {U) 
might be responsible not only for the formation of 
the HF state but also for the occurrence of HF su- 
perconductivity. 

2. In the superconducting state, for CeCu2Si2, 
CeCoIns and UPd2Al3, the spin part in the Knight 
shift deduced from or Kt^ is in excellent agree- 
ment with the reduction in the observed Knight 
shifts 6K°'"^ below T^, giving unambiguous evi- 
dence that the heavy quasi-particles form the spin- 
singlet Cooper pairs. By contrast, based on the 
present analyses, UPts and UNi2Al3 are concluded 
to belong to a class of odd-parity (spin-triplet) pair- 
ing state. 



The present NMR analysis has shed light on the semi- 
quantitative estimate for the quasi-particles susceptibil- 
ity in the HF systems. Especially from the present anal- 
ysis based on the Fermi-liquid theory, it is suggested that 
the renormalizcd HF ground state has a Kramers degen- 
eracy through the c — / mixing even for U-based HFS's. 
However, we are still at a long way from a microscopic 
imderstanding how to represent the heavy-Fermi liquid 
state enhanced through the hybridization with a singlet 
CEF ground state. It is needed to be addressed by fur- 
ther experimental and theoretical efforts in future. 
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